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ABSTRACT 



o 

^ ' We present estimates of the basic physical properties (size and albedo) of (90377) Sedna, a prominent member of the 

detached trans-Neptunian object population and the recently discovered scattered disk object 2010 EK139, based on the 
recent observations acquired with the Herschel Space Observatory*, within the "TNOs are Cool!" key programme. Our 
modeling of the thermal measurements shows that both objects have larger albedos and smaller sizes than the previous 
expectations, thus their surfaces might be covered by ices in a significantly larger fraction. The derived diameter of 
Sedna and 2010 EK139 are 995 ± 80 km and 470^ 3 q km, while the respective geometric albedos are py = 0.32 ± 0.06 and 
0.25^q'o5. These estimates are based on thermophysical model techniques. 

Key words. Kuiper belt objects: (90377) Sedna, 2010 EKi3g - Radiation mechanisms: thermal - Techniques: photometric 

o 

, 1. Introduction Sedna is a prominent member of the detached objects, 

■ that is often classified as an inner Oort-cloud object. Until 

The Herschel Space Observatory (jPilbratt et all l2010l ) al- now, no accurate measurements of the diameter and albedo 
lows the detection of thermal radiation from several trans- have been ava ilable for this object. Both direct imaging 
Neptunian objects (TNOs) at the precision level of < (iBrownl . I2008D and upper limits to the thermal radiation 
lmJy. Since the expected fluxes around the peak of the using the Spitzer Spac e Telescope (2.4 mJy at 70 ^m, see 
H . spectral energy distribution (SED) significantly exceed this IStansberrv et all . 120081) have yielded an upper limit of re 
precision, Herschel provides a great opportunity to char- 1670km for its diameter (within 97% confidence), 
acterize TNOs and obtain basic thermophysical informa- 2010EK 139 has been discovered in 2010 by 
tion. In this work, we present recent measurements of the ISheppard et all ([20lTh in the course of a southern 
prominent objects (90377) Sedna and 2010 EKi 39 using the Galactic plane survey. Prediscovery observations date back 
Photodetector Array Camera and Spectrometer instrument to 2 002, allowing for a relatively accurate orbit determi- 
(PACSl Poglitsch et aU |2fll0) on board the Herschel Space nation This places 2 010EK 139 among the scattered disk 
Observatory. These observations are part of the "TNOs are ob j e cts. 2010EK 139 orbits the Sun on an eccentric orbit 
Cool!: a survey of the trans-Neptunian region" Open Time ( e ~ 53) and has a perihelion distance of q re 32.5 AU. 
Key Program [Mu^ et aL l 12Q0JI 120101 | Lellouch et al. |, In addit i on7 2010EK 139 is a suspected member of the 2:7 
[20101; | Lim et al. |, | 2010D . resonance groupQ. We note that a more complete sample 

of SDOs/dctac hed objects observed wit h Herschcl/PACS 
is presented bv lSantos-Sanz et all (12012ft . 



X 



* Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consor- 

tia and with important participation from NASA. 1 http://boulder.swri.edu/~buie/kbo/astrom/10EK139.html 
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Table 1. Summary of Herschel observations of Sedna and 
2OIOEK139. The columns are: i) observation identifier, ii) 
date and time, iii) scan angle direction with respect to the 
detector array, and iv) filter configuration. 



70 fim 



100 /mi 



160 (im 



OBSID 


Date & Time (UT) 


Angle 


Filter 


Sedna 


1342202227 


2010-08-06 10:55:17 


70° 


B/R 


1342202228 


2010-08-06 11:19:54 


110° 


B/R 


1342202229 


2010-08-06 11:44:31 


70° 


G/R 


1342202230 


2010-08-06 12:09:08 


110° 


G/R 


1342202306 


2010-08-09 08:11:37 


70° 


B/R 


1342202307 


2010-08-09 08:36:14 


110° 


B/R 


1342202308 


2010-08-09 09:00:51 


70° 


G/R 


1342202309 


2010-08-09 09:25:28 


110° 


G/R 


2010 EK139 


1342211418 


2010-12-23 07:04:30 


70° 


B/R 


1342211419 


2010-12-23 07:15:01 


110° 


B/R 


1342211420 


2010-12-23 07:25:32 


70° 


G/R 


1342211421 


2010-12-23 07:36:03 


110° 


G/R 


1342211524 


2010-12-23 19:58:27 


70° 


B/R 


1342211525 


2010-12-23 20:08:58 


110° 


B/R 


1342211526 


2010-12-23 20:19:29 


70° 


G/R 


1342211527 


2010-12-23 20:30:00 


110° 


G/R 



2. Observations, data reduction and photometry 

Sedna was observed by Herschel/PACS in two visits: the 
first started on 2010 August 6, 10:55:17 UTC and a follow- 
up started on 2010 August 9, 08:11:37 UTC, both taking 
place during the Routine Science Phase obs ervation series 
of th e "TNO's are Cool!" key programme (jMuller et all . 
12009ft . 2OIOEK139 was also observed by Herschel/PACS in 
two visits, the first started on 2010 December 23, 07:04:30 
UTC, and a follow-up started the same day, 19:58:27 UTC. 
Herschel/PACS observed Sedna and 2010EKi 39 for « 3.14 
and « 1.26 hours, respectively. For both objects, we used 
both the blue/red (70/160 fjm) and green/red (100/160 /xm) 
channel combinations. The actual details of these observa- 
tions are summarized in Table [TJ 

Raw observational data were reduced using the Herschel 
Interactive Processing Environment (HIPE0, see also Ott, 
2010) an d the processing scripts are similar to the ones em - 
pl oyed inlMommert et al](|2012[ ). ISantos-Sanz et all (|2012t) . 
or lVilenius et al.l ( 2012[ ). For each observation, these scripts 
create a pair of maps, one for the blue or green chan- 
nel and one for the red channel. The maps have an effec- 
tive pixel size of l'.'l, l'/4, and 2'.'1, for the blue, green, 
and red filters, respectively: these pixel sizes are set to 
sample the respective point spread functions (PSFs) prop- 
erly. Data frames were selected by the actual scan speed 
(10"/sec < speed < 25"/sec) of the spacecraft, which max- 
imized the effective usage of the detector and yielded signifi- 
cantly higher signal-to- noise (S/N) ratios than the standard 
setting (approximately 20"/sec). 

Since the apparent motion of Sedna and 2OIOEK139 be- 
tween the two visits (15 — 35 map pixels, depending on 
the actual filter) is relatively large compared to the PSF 
but small compared to the detector size, the location of 



2 Data presented in this paper were analyzed using "HIPE", 
a joint development by the Herschel Science Ground Segment 
Consortium, consisting of ESA, the NASA Herschel Science 
Ce nter, and the HIFI, P ACS an d SPIRE consortia mem bers, 
see |http://herschel.esac.esa.int/DpHipeContributors.shtml| 
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Fig. 1. Image stamps showing the combined maps of Sedna 
(upper panels) and 2010EKi3g (lower panels) in the 70 /mi 
(blue), 100 /im (green), and 160 /im (red) channels. Each 
stamp covers an area of 56" x 56", while the tick marks on 
the axes show the relative positions in pixels. The effective 
beam size (i.e. the circle with a diameter corresponding to 
the full width at half magnitude) is also displayed in the 
lower-left corners of the stamps. 

Table 2. Thermal fluxes of Sedna and 2010EKi 39 derived 
from our Herschel measurements. 



Object 


Band 


A 


Flux 


Sedna 


B 


70 /im 


1.8±0.7mJy 




G 


100 /mi 


4.2±0.9mJy 




R 


160 /im 


2.7 ± 1.3 mJy 


2010EK 139 


B 


70 /im 


17.4 ± 1.1 mJy 




G 


100 /im 


16.3 ± 1.4 mJy 




R 


160 /im 


11.9 ± 1.8 mJy 



the target in the first visit can simply be used as a back- 
ground area on the maps of the second visit and vice versa. 
Owing to t he satellite pointing u ncertainty that is about a 
few arcsec (|Poglitsch et aUl2010j) . we derived the true map- 
center displacements using the red channel maps - on which 
the background confusion is the strongest - as follows. By 
varying the proper motion vector between the two visits, 
we computed the cross-correlation residuals for each trial 
vector. By minimizing the residuals, we obtained a more 
precise value for the shift between the visits and the pho- 
tometry of combined maps was found to be more reliable. 
Since simple averaging the registered maps does not cancel 
the background confusion noise, we em ployed background 
remov al t echniques as it is descr ibed in Sa ntos-Sanz et all 
(|2012l ) or iMommert et al.l ()2012h . The maps on which the 
photometry was then performed are shown in Fig. [TJ 

Regardless of the background structures, in the sub- 
tracted and combined maps the only expected source is 
the TNO itself and this source can be treated as an iso- 
lated point source. We estimated the fluxes and their un- 
certainties using (1) a single aperture that maximizes the 
expected S/N ratio; (2) the aperture growth curve method 
and imp lanted artificial sources in a Monte Carlo fashion 
(see e.g. iSantos-Sanz et al.L l2012h ; and (3) we also checked 
the individual (non-combined) maps on which they had suf- 
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Fig. 2. Spectral energy distribution of Sedna (left) and 2OIOEK139 (right) in the far-infrared region, based on 
Herschel/PACS measurements. Superimposed are the best-fit TPM (solid lines) and STM curves with floating beaming 
parameter (dashed lines). 



Table 3. Orbital and optical data for Sedna and 
2OIOEK139 at the time of the Herschel observations. The 
parameters r and A denote the heliocentric distance and 
the distance from Herschel, a is the phase angle, and Hy is 
the absolute visual magnitude, which is available from the 
literature. 



Object 


Quantity 


Value 


Sedna 


r 


87.43 AU 




A 


87.56 AU 




a 


0?7 




H v 


+1.83 ±0.05 


2OIOEK139 


r 


39.08 AU 




A 


39.50 AU 







1?3 






+3.80 ±0.10 



fic ient S/N ratio. For a m ore detailed description , we re fer 
to lMommert et al.l (1201 2D and ISantos-Sanz et al.l (|2012l ). 

We found that all three methods yielded the same fluxes 
and uncertainties for each channel. The individual analy- 
sis of maps for 2OIOEK139 also showed consistent results. 
Therefore, we accepted the means of all measurements per 
object and filter as final fluxes (see Table [5]) and used them 
for thermal mod eling. We note here th at the color correc- 
tions provided in lPoglitsch et all (|2010l ) are negligible: it is 
nearly or less than 1 per cent for 2010 EKi3g and for Sedna, 
it is 6 per cent in the 70 /jm channel and less than 3 per cent 
in the other longer wavelength channels, so almost less by 
a magnitude than the relative photometric uncertainties in 
all of these cases. 



3. Thermal properties 

The basic physical properties of Sedna and 2OIOEK139 
were e stimated by a hybri d standard thermal mo del 
(STM, iLebofskv et all . 119861: IStansberrv emi . I2008D in 
which the beaming parameter is adjustable and the 

aster oid th ermop hysical model (TPM, see iLagerrosl 

119961 I1997L 119981 ). The absolute magnitudes of the re- 
flected sunlig ht from these TNOs are available from 
the literature (Rabinowitz. Schaefer fc Tourtellottd . l2007t 



ISheppard et all 120111 ). In addition for 2010EKi- ; q. we 



conservatively increased the formal uncertainty (0.03 mag, 
based on MPC data) up to 0.10 mag: we took into account 
the possible omission of the phase angle corrections and 
also added quadratically an av erage TNO lightcury e am- 
plitude of 0.088 mag (based on iDuffard et all . 120091 ). The 
employed geometric parameters and absolute magnitudes 
are summarized in Table [3] 

The hybrid STM predicts thermal fluxes from the geo- 
metric albedo, diameter, and beaming parameter and these 
fluxes can be computed for arbitrary solar and geocen- 
tric distances. Hybrid STM provides reliable estimates only 
for small phase angles (via a simple form of phase angle 
corrections), although owing to the distances of these ob- 
jects, this estimate is fairly sufficient in our cases. To es- 
timate the physical parameters and their respective uncer- 
tainties, we used a Monte Carlo (MC) approach by varying 
the fit input around their mean with the standard devi- 
ation equal to their respective uncertainty. For both tar- 
gets, we used the fixed- 77 approach for the beaming pa- 
rameter, i.e. taking 77 = 1.14 ± 0.15 for each MC step. 
This mean value of and scatte r in th e beaming factor are 
taken from ISantos-Sanz et all (|2012| ) and seem to be an 
acceptable approach for TNOs. To estimate the phase in- 
tegral q, i.e. the ratio of the Bond to geometric albedo 
(i.e. A = qpv), we employed an iterative approach. First, 
the phase integral is computed for unity slope parameter 
(G = 1 , i.e. q = 0.2 9 + 0.6 8G), and then refined using 
cq. 1 of iBrucker et al.l (|2009l ) until convergence. This pro- 
cedure applied for hybrid STM yielded the diameter, geo- 
metric albedo, and slope parameter of D = 1060 ± 100 km, 
pv = 0.290 ± 0.061, and G = 0.42 ± 0.04 for Sedna and 
D = 535 ±30 km, p Y = 0.187+0.027, and G = 0.37 ±0.03 
for 2OIOEK139, respectively. We repeated the similar pro- 
cedure by allowing the beaming parameter 77 to vary. This 
analysis yielded r\ = 0.95 ± 0.43 for Sedna, with the cor- 
responding diameter and albedo of D = 990 ± 95 km and 
pv = 0.336 ± 0.072. For 2010EKi 39 , the best-fit value of 
the beaming parameter is somewhat smaller, 0.70 ± 0.31, 
while the diameter and albedo values arc D = 450 ± 35 km 
and p v = 0.261 ± 0.047. 

In the case of Sedna, we note that the 

l inear phase coefficient (3 = 0.151 ± 0.033 

(jRabinowitz. Schaefer fe Tourtellottel . l2007t) would imply 
a phase integral of q = 0.891^' 29, assuming the same phase 
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behavior over the whole phase angle range. Although the 
phase curve is known for very small domains (a < - 6° , se e 
also fig. 2 of iRabinowitz. Schaefer k Tourtellottel l2007ft . 
this value broadly agrees, within a nearly 1-a difference 
from the phase integral of q = 0.59 ± 0.03 as implied by 
the radiometric albedo and the Brucker formula. 

The results of the TPM estimates were the following. 
For Sedna, we used the rotation period of w f 0.27 hours 
(|Gaudi et all 120051 ) and assumed an equator-on rotation 
and the most favorable solution for the thermal inertia 
was found to be 0.2 J m~ 2 K _1 s -1 / 2 , which corresponds 
to the diameter of D = 995 ± 80 km and the geometric 
albedo of p Y = 0.32 ± 0.06. For 2010EKi 39 , we assumed 
a period of 12 hours and equator-on rotation and found 
that this object also requires a very low thermal inertia, 
0.1Jm _2 K _1 s -1 / 2 . This may change slightly if the ro- 
tation period and the spin vector orientation were very 
different, although all feasible solutions put the thermal 
inertia below f .0 J m~ 2 K _1 s -1 / 2 . Our best fit yielded a 
diameter of D = 470+?q km and a geometric albedo of 
Pv = 0.251qq 2 , which do not differ significantly from the 
hybrid STM model results in which the beaming param- 
eter was also varied. We note that here the uncertainties 
include both the statistical errors and the amb i guitie s in 
the rotation parameters (see also iMiiller et all 201 ll for 
more details about this modeling). 

In Fig. [2j we displayed the far-infrared SEDs for these 
two objects as it is estimated from the hybrid STM and 
TPM fitting and our best-fit data and the floating rj val- 
ue s. We also no t e tha t the rotation period of Sedna found 
by iGaudi et al.l (|2005ft corresponds to a peak-to-pcak am- 
plitude of 0.02 mag. The small amplitude does not change 
the reliability of the thermal modeling and the correspond- 
ing shape effects are not relevant to the size determination. 

4. Discussion 

We have estimated the sizes and surface albedos for 
the trans- Neptunian objects Sedna and 2OIOEK139 us- 
ing recent observations of their thermal emission at 
70/100/160 nm with Herschel/PACS. On the basis of ear- 
lier Spitzer measurements, S edna had already only a n up- 
per limit to its size estimate (|Stansberrv et all 1 2 00 8ft . Our 
analysis has shown that for Sedna and 2OIOEK139, the re- 
spective geometric albedos are p\r = 0.32 ± 0.06 and py = 



0.25; 



-0.02 

-0.05 



thus both objects have brighter surfaces t han 



the average TNO population (jStansberrv et all 20081) or 



SDO s /detached population (table 5 in ISantos-Sanz et all 
12012ft . We note that the albedos of Sedna and 2010EKi 39 
closely match th ose of detached objects in fig. 4a (p\r 
vs. diameter) in Santos-Sanz et al. (12012ft. According to 



I Schaller k Brown ( 2007ft and iBrown. Burgasser k Fraserl 
(|201lft . even with these newly derived parameters, Sedna 
lies in the region in which volatiles are expected to be re- 
tained in the surface (see Fig. 1 in that paper for an equiv- 
alent temperature of 20 ±2K), hence one can also expect a 
brigh ter surface (see also lBarucci et all 120051 ; lEmerv et all 
12007ft . Sedna is currently approaching its perihelion. Thus, 
if the brightness of the surface were changing owing to 
the ongoing sublimation of ices, it might be detectable in 
the variation in the absolute magnitude on a timescale of 
decades. In contrast, 2010 EK139 falls in the region in which 
volatiles should have been lost (using 48 ± 3 K for equiv- 
alent temperature). However, objects of this size can have 



such a h igh albedo if water ice is present on the surface 
(see e.g. Barkume et all 120061 : iRagozzine k Brownl . 120091 : 
iDumas et all 2011ft . The presence of water could be tested 
by measuring the intrinsic color to see whether it is bluish 
(iBrownl . [2008J). As they lack known satellites, we do not 
know the masses and hence the surface gravity and es- 
cape velocities of these objects that could place tighter con- 
straints on the surface properties. 
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